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Abstract. Let A* > denote the largest possible value of A such that 

A^u = A(l + u)P in 



CN ■ \ u = ^ = on dM 

' has a solution, where B is the unit ball in i?" centered at the origin, 

Q> . p > and n is the exterior unit normal vector. We show that for 

^ I A = A* this problem possesses a unique weak solution u*, called the 

extremal solution. We prove that u* is singular when n > 13 for p large 

4 

enough, in which case u*{x) < r p-i — 1 on the unit ball and actually 
solve part of the open problem which ^ left. 

^ . 1. Introduction and results 

I In the previous two decades, positive solutions to the second order semilinear elliptic 

problem 

-Au = Xg{u) in f2, 

n = on n, ^ ' 

have attracted a lot of interest, see e.g. [IHS] and references therein. Here, we only mention 
the work by Joseph and Lundgren [2]. In their well known work, Joseph and Lundgren 
gave a complete characterization of all positive solutions of (1.1) in the case g{u) = 
or g{u) = (1 + auy, ap > 0, X > and Q is unit ball in i?". In particular, they found a 
remarkable phenomenon for g{u) = e" and n > 2: either (1.1) has at most one solution 
for each A or there is a value of A for which infinitely many solutions exist. In the case of a 
power nonlinearity the same alternative is valid if > 3 and p ^ (1, ^^]. As a subsequent 
step, P.L. Lions ( [3], section 4.2 (c)) suggests to study positive solutions to systems of 
semilinear elliptic equations. So it is an important task to gain a deeper understanding 
for related higher order problems. 
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In this paper we study a semilinear equation involving the bilaplacian operator and a 
power type nonhnearity 

A^u = A(l + u)P in B, 



u = ^ = on dM 

an 



where B C i?" is the unit ball, A > is an eigenvalue parameter, n > 5 and p > ^^l- The 
subcritical case p < is by now "folklore" , where existence and multiplicity results are 
easily established by means of variational methods. For the critical case p = (under 
Navier boundary conditions), we refer to [6]. Recently, a lot of research on supercritical 
case , i.e., p > has been done and many beautiful important results have been proved. 
In what follows, we will summarize some of the results obtained by [UIH] . For convenience, 
we introduce the following notions: 

Definition 1.1. We say that u G Lp{B) is a solution of (1.2) if u > and if for all 

if e C^(B) with ifldM = \V(p\\dM = one has 



/ uA^ifdx = A / (1 + ufipdx. 
Jb Jb 



We call u singular if u ^ L°°(B), and regular if u G L'^ 

4 

A radial singular solution u = u{r) of (1.2) is called weakly singular if limr_j.o ^^-^^(r) G 
[0, oo] exists. 

Note that by standard regularity theory for the biharmonic operator, any regular 
solution u of (1.2) satisfies u G C°°(B). Note also that by the positivity preserving 
property of in the ball (see [TU]) any solution of (1.2) is positive, see also [TT] for a 
generalized statement. This property is known to fail in general domains. For this reason, 
we restrict ourselves to ball. Hence, the sub-and super-solution method applies as well as 
monotone iterative procedures. 

Definition 1.2. We call a solution u of (1-2) minimal if u < v a.e. in B for any further 
solution V of (1.2) 

We also denote by Ai > the first eigenvalue for the biharmonic operator with Dirichlet 
boundary conditions 

A'^ipi = Aiv^i in B, 
^^ = 9^ = on dM. 

It is known from the positivity preserving property and Jentzsch's (or Krein-Rutman's) 
theorem that Ai is isolated and simple and the corresponding eigenfunction ipi don't 
change sign. 

Definition 1.3. We say a weak solution of (1.2) is stable (resp. semi-stable) if 



inf{ / {A^f-pX f + : G H',{M), = 1} 



is positive (resp. non-negative). 

To illuminate the motivations of this paper in detail, we need the following notations 
which will be used throughout the paper. Set 

8(j9 + i) r 2(p + i)i r 4p 

Ko = — n — n , 

p — 1 p — 1 p — 1 



and 

n + 2 - ^/A + ^^^^IV^T^ , 
Pc = tor n > 6, 

with Hn = {n{n — 4)/4)^ and the number pc such that when p = Pc then 

4 4 4 4 

+ 4)( + 2){n - 2 -){n - 4 -) = H^. 



p — 1 p — 1 p — 1 p — 1 

Now we summarize some of the well-known results as follows 



Kq, — ) such that: 



Theorem A There exists A* G 

(i) For X G (0, A*), (1-2) admits a minimal stable regular solution, denoted by ux. This 
solution is radially symmetric and strictly decreasing in r = \x\. 

(a) For A = A*, (1-2) admits at least one not necessarily bounded solution, which is called 
extremal solution u* 

(Hi) For A > A*, (1-2) admits no (not even singular ) solutions 
Theorem B /P./. Assume that 

< p < Pc ij n > 16, <j9<oo ij b < n < 12 

n — A n — A 

Then, u* is regular. 

From Theorem B, we know that the extremal solution of (1.2) is regular for a certain 
range of p and n. At the same time, they left a open problem: if 

n> 13 and p > Pc, 

is u* singular? 

In this paper, by constructing a semi-stable singular Hq(E)— weak sub-solution of 
(1.2), we prove that, if p is large enough, the extremal solution is singular for dimensions 
n > 13 and complete part of the above open problem. Our result is stated as 
follows: 

Theorem 1.1. There exists Pq > 1 large enough such that for p > Pq, the unique extremal 

4 

solution of (1.2) is singular for dimensions n > 13, in which case u* < \x\ p-^ —1 on the 
unit ball. 



From the technical point of view, one of the obstacle is the well-known difficulty of 
extracting energy estimates for solutions of fourth order problems from their stability 
properties. Besides, for the corresponding second order problem (1.1), the starting point 
was an explicit singular solution for a suitable eigenvalue parameter A which turned out to 
play a fundamental role for the shape of the corresponding bifurcation diagram, see [12j . 
When turning to the biharmonic problem (1.2) the second boundary condition |^ = 
prevents to find an explicit singular solution. This means that the method used to analyze 
the regularity of the extremal solution for second order problem could not carry to the 
corresponding problem for (1.2). In this paper, we, in order to overcome the second 
obstacle, use improved and non standard Hardy-Rellich inequalities recently established 
by Ghoussoub-Moradifam in [13] to construct a semi-stable singular Hq(M)— weak sub- 
solution of (1.2). 



This paper is organized as follows. In the next section, some preliminaries are reviewed. 
In Section 3, we will show that the extremal solution u* in dimensions > 13 is singular 
by constructing a semi-stable singular Hq{M)— weak sub-solution of (1.2). 

2. Preliminaries 

First we give some comparison principles which will be used throughout the paper. 
Lemma 2.1. (Boggio's principle, 110]) If u E C^(]B/j) satisfies 

( A^u > in Mr, 
\u = ^ = on dMn, 

then u > in Mr. 

Lemma 2.2. Let u G L}{Mr) and suppose that 

I mAV > 

for all G C'*(B/j) such that (f > in Mr, fldMn = |f laB^ = 0. Then u > in Mr. 
Moreover u = or u > a.e., in Mr. 

For a proof see Lemma 17 in [11]. 

Lemma 2.3. If u & H'^{Mr) is radial, A'^u > in Mr in the weak sense, that is 

/ AuA<^ > G C^{Mr), > 

and wlaBfl > 0, ^Iomr < 0, then u>0 in Mr. 

Proof. For the sake of completeness, we include a brief proof here. We only deal with 
the case R = 1 for simplicity. Solve 

f A^Mi = A^M in B 
I m = |ii = on dM 

in the sense Ui G H^{M) and J^AuiAif = J^AuAif for all (f G C^(B). Then mi > in 
B by lemma 2.2. 

Let U2 = u — Ui so that A^M2 = in B. Define / = Au2. Then A/ = in B and since 
/ is radial we find that / is a constant. It follows that U2 = ar"^ + b. Using the boundary 
conditions we deduce a -|- > and a < 0, which imply U2 >0. 

Now we give a notion of Hq{M)- weak solutions, which is an intermediate class between 
classical and weak solutions. 

Definition 2.1. We say that u is a H^(M)- weak solution of (1.2) if (1 + uY G ^^(B) 
and if 

I AuAcj) = A /" 0(1 + uf, V0 G C\M) n Hl{M). 
Jb Jb 

We say that u is a Hq(M)- weak super- solution (resp. Hq(M)- weak sub-solution) of (1.2) 
if for (p > the equality is replaced with > (resp.<) and u > (resp. <), |^ < (resp. 
>) on dM. 

We also need the following comparison principle. 



Lemma 2.4. Letui,U2 G H'^{^r) with (1 + Ui)^, (1 + ^2)^ G L^(Bij). Assume that Ui is 
stable and 

A^Mi < A(l + m Mr 
in the i7^(B/j)— weak sense, i.e., 

I AuiA^<xf y4>eC^(Mn),4>>0 (2.1) 

anc? A^U2 > A(l + ^2)^ m B^j in the similar weak sense. Suppose also 

dui du2 1 

Then 

Ui < U2 in B/j. 

Proof. Define u := Ui — U2- Then by the Moreau decomposition [T3] for the biharmonic 

operator, there exist uji,uj2 G HliMji), with u = uji + (jJ2,^i > a.e., A^u)2 < in the 
HI{Mb)— weak sense and 

/ AwiAw2 = 0. 

By Lemma 1.1, we have that a;2 < a.e. in B^j. 
Given now < G C^(Bij), we have that 

/ Au;A<^ < A / (/(Ml) - /(n2))<^, 

where f{u) = (1 + uY- Since u is semi-stable and by density one has 

A / f\u)ul < A / (Aa;i)2 = A / Ac^Awi < A / (/(m) - /(n2))a;i. 

JBfl JKr JEr JEr 

Since cji > one also has 

/ nu)uu, < j (/(Ml) - f{u2))uJi 

JMr JEr 
which once re-arrange gives 

/ /c^i>0, 
Jer 

where f{ui) = f{ui) — f{u2) — f'{ui){ui — U2). The strict convexity of / gives / < and 
/ < whenever u ^ U . Since cji > a.e. in B^j, one sees that u; < a.e. in B/j. The 
inequahty ui < U2 a.e. in B/j is then estabhshed. 

The folfowing variant of lemma 2.4 also holds: 

Lemma 2.5. Let Ui,U2 G H^{Mr) be radial with (1 + UiY, (1 + ^2)^ G L^(Br). Assume 
A^ui < A(l + Mi)P m Mfi in the sense of (2.1) and and A?U2 > A(l + ^2)^ in B^j. 
Suppose uild^R < U2\dMR and ^Isbr > ^{omr O'^d suppose also that Ui is semi-stable. 
Then Ui < U2 in B^j. 



Proof. We solve for u e H^{M) such that 

f AuA(f)^ I A{ui-U2)A(f) V(/) e Co°°(Br). 

By Lemma 2.3 it follows that u > ui — U2. Next we apply the Moreau decomposition to 
M, that isix = w + v with w,v e H^{Mr),w > 0, A^v < in B^j and AwAv = 0. Then 
the argument follows that of Lemma 2.4. 

Lemma 2.6. Let u he a semi-stable Hq{M)— weak solution of (1.2). Assume U is a 
Hq{E)— super-solution of (1.2). Then if u is a classical solution and Hi{u) = 0, we have 
u = U. 

Proof. Since m is a classical solution, it is easy to see that the infimum in fj,i{u) is attained 
at some (p. The function (p is then the first eigenfunction of — Xf'{u) in Hq(B), where 
f{u) = (1 + uy. Now we show that is of fixed sign. Using the Moreau decomposition, 
one has = 0i + 02 where 0^ G H^iM) for i = 1, 2, 01 > 0, A0iA02 = and A^^a < 
in the Hq{M)— weak sense. If changes sign, then 0i ^ and 02 < in B. We can write 
now: 

^ , , ^ /b(A(0i - 02))^ - Af (^)(0i - 02)^ ^ ^(A0)^ - Xf{u)<p' 

in view of 0i02 < — 0i02 in a set of positive measure, leading to a contradiction. 

So we can assume > 0, and by the Boggio's principle we have > in B. For 
< t < 1 define 



g{t) = / A{tU + (1 - t)M)A0 - A / f{tU + (1 - t)M)0, 
Jb Jm 

where is the above first eigenfunction. Since / is convex one sees that 

9{t) >X [ [tf{U) + (1 - t)f{u) - f{tU + (1 - t)u)]<P > 

for every t > 0. Since g{0) — and 

g'{Q)= I A{U -u)A(j)- Xf'{u){U -u)(p = Q 
Jb 

we get that 

/(0) = -A / nu){U-uf<p>0. 
Jb 

Since f"{u)(f) > in B, we finally get that U — u a.e. in B. 
Prom this lemma, we immediately obtain: 

Corollary 2.1 (i) When u* is a classical solution, then Hi{u*) — and u* is the unique 
Hl{B>)— weak solution of (1.2); 

(a) If V is a singular semi-stable Hq{M)— weak solution of (1.2), then v — u* and A = A*. 

Proof, (i) Since the function u* is a classical solution, and by the Implicit Function 
Theorem we have that ^J^l{u*) = to prevent the continuation of the minimal branch 
beyond A*. By Lemma 2.4, u* is then the unique i?Q(B)— weak solution of (1.2). 



(ii) Assume now that f is a singular semi-stable Hq{M)— weak solution of (1-2). If 
A < A*, then by the uniqueness of the semi-stable solution, we have v = u\. So v is not 
singular and a contradiction arises. By Theorem A (iii) we have that A = A*. Since v is 
a semi-stable Hq{M)— weak solution of (1.2) and u* is a Hq{W)— weak super-solution of 
(1.2), we can apply Lemma 2.4 to get v < u* a.e. in B. Since u* is a semi-stable solution 
too, we can reverse the roles of v and u* in Lemma 2.5 to see that v > u* a.e. in B. So 
equality v = u* holds and the proof is complete. □ 

3 Proof of Theorem 1.1 

Inspired by the work of [16j, we will first show the following upper bound on u* 

4 

Lemma 3.1. If n > 13 and p > Pc, then u* < \x\ p-'^ — 1. 



4 



Proof. Recall from Theorem A that Kq < A*. We now claim that u\ < u := \x\ p-i 
for all A G {Kq, A*). Indeed, fix such a A and assume by contradiction that 

Ri := inf{0 < R < 1 : Ux < u in the interval {R, 1)} > 0. 

From the boundary conditions, one has that 

Ux{r) < u{r) as r — 7- 1~. 

Hence, 

< i?i < l,uxiRi) = uiRi) and u'x{Ri) < u'{Ri). 
Now consider the following problem 

A^u = Ko{l + uy in B^^,; 



u = ux{Ri) on dMji,; 

|^ = <(i?i) on dMn,. 



Then ux is a super-solution to above problem while m is a sub-solution to the same 
problem. Moreover for n > 13, we have 

pKo<H^:= ^ ' 



16 

and 



/ (A0)2 >Hr, f ^rfx > pKq I (1 + 



So u is semi-stable and we deduce that ux > u hj the Lemma 2.4, and a contradiction 
arises in view of the fact 

\u\\l^(Mr,) < oo, and \u\l^(Br,) = oo. 

The proof is done. □ 

In order to prove Theorem 1.1, we will need a suitable Hardy- Rellich type inequality 
which was established by Ghoussoub-Moradifam in [13]. It is stated as follows: 

Lemma 3.2. Let n > 5 and B be the unit ball in R". Then there exists C > 0, such that 
the following improved Hardy-Rellich inequality holds for all ip G Hq{ 

I {AifYdx > / ^rfx + C [ ^Hx. 

Jm 16 Jb fI Jm 



Lemma 3.3. Let n > 5 and B be the unit ball in M". Then the following improved 
Hardy-Rellich inequality holds for all ip G H'q{M): 



I 

Jb. 



(n - 2)2(n - 4)2 f ip'^dx 



x2, (n-2)\n-4:f f 

(Aip)'dx > ^ L / 

16 Jb 



(|a;|2 -0.9|a;|t+i)(|a;|2 - \x\^) 



. (!i^iK!i^ / (3.0) 



As a consequence, the following improvement of the classical Hardy-Rellich inequality 
holds: 

We now give the following lemma which is crucial for the proof of the Theorem 1.1. 

Lemma 3.4. Suppose there exist A' > and a radial function u e H^{M)r)W^^^{B\{0}) 
such that u L°°(B) and 

A\ < A'(l + uf for < r < 1; u{l) = u'{l) = 

and 

pf3 [ ip\l + uy+^ < [ {Aiff for all if e H^{M) 
Jb Jb 

for either ^ > A' or ^ = A' = Then u* is singular and 

A* < A'. (3.1) 

Proof. First, we prove A* < A'. Note that the stability and u e (B \ {0}) yield to 
(1 + uY e -£''^(B), we easily see that u is a, weak sub-solution of (1.2). If now A' < A*, by 
Lemma 2.5, u would necessarily be below the minimal solution uy, which is a contraction 
since u is singular while uy is regular. 

Suppose first that (3 = X' — ^ and that n > 13. From the above we have A* < -y^, wc 
get from Lemma 3.1 and the improved Hardy-Rellich inequality that there exists o > 



so that for all e 



I (A0)2 - pX* U\l + uy^' > [ (A0)2 -H^ [ 
Jb Jb Jb Jb 



/>2 

> c 



X 



4 



It follows that IjLi{u*) > and u* must therefore be singular since otherwise, one could 
use the Implicit Function Theorem to continue the minimal branch beyond A* 

Suppose now that /3 > X' and let ^ < 7 < 1 and a := (^)^ and define u :— 
a~^{l + u) — 1. We claim that 

u* <u in B. (3.2) 
To prove this, we shall show that for every A < A* 

ux<u in B. (3.3) 



Indeed, we have 

A^u = aA\ < aX'{l + uf = aP+^X'{l + uf. 

Now by the choice of a, we have a'^'^^X' < X*. To prove (3.3), it suffices to prove it for 
aP+^X' < A < A*. Fix such A and assume that (3.3) is not true. Then 

A = {0 < < l\uxiR) > 

is non-empty. Since u{l) = — 1 > = ma(1), we have < it!i < l,ux{Ri) — u{Ri), 
and u\{Ri) < u'{Ri). Now consider the following problem 

A^u^X{l + u)P in Mr^, 
u = ux{Ri) on 9Biji, 

g = <(i?i) on OMr,. 

Then ux is a solution to above problem while u is a. sub-solution to the same problem. 
Moreover u is stable since A < A* and 

pA(l + uf+' < pX*a-^^^'\l + uf^' = pX'^-\l + uy+' < pPi{1 + uy+\ 

we deduce u < uxin Mpt-^ , which is impossible, since u is singular while ux is smooth. This 
establishes (3.2). From (3.2) and the above inequalities, we have 



Thus 



pA*(l + uy+^ < pA'7"'(l + uf+^ < pA(1 + uy+\ 
/„(A^)2-pAV(l + «T+' 



inf 

<pecg°(u) 



> 0. 



This is not possible if u* is a smooth function by the Imphcit Theorem. 



□ 



Proof Theorem 1.1 Uniqueness and the upper bound estimate of the extremal solution 
u* have been proven by Corollary 3.1 and Lemma 3.1. Now we only prove that u* is a 
singular solution of (1.1) for n > 13, in order to achieve this, we shall find a singular 
if— weak sub-solution of (1.1), denote by ujm{r), which is stable, according to the Lemma 
3.4. 



Choosing 



Um = air f-i + a2r 



since uj{l) = uj'{l) = 0, we have 



ai 



1, i^o = 



S{p + 1) 



P 



n — 



2{p+l) 



P 



n — 



4p 



P 



III 



a2 



4 
p-1 



p-1 



p-1 



For any m fixed, when p +oo, we have 

4 



Oi 



1 - 



(p — l)m 



+ 0{p ^), 02 = 1 - Ol 



{p — l)m 



+ o{p-^) 



and 



^^^ 8(»-2)(,.-4) ^ 

P 



Note that 



4p 



m— 4 



4p 



A'ii'o(air + a2r"^Y - aiKor p+i - a2Kir 



m— 4 



4p 



Ap 



A'(ai + aar'^+F^)^' - ai - asi^ii^o"'^^^"""^] 
A'(ai + a2r"^+^)^ - ai - a2i^iii:o"V"^-^] 



4p 



(3.4) 



with 



H{x) = (oi + 02^)^ [oi + a2i^ii^o'^^] > -^i = ^("^ - 2)(m + n - 2)(m + n - 4) (3.5) 
(1) Let m — 2 and n > 32, then we can prove that 



sup if (x) = if(0) = 

[0,1] 



e as p 



+00. 



So (3.4) > is vahd as long as 



X — e . 



At the same time, we have ( since ai + a2r p-i < Oi + 02 < 1 in [0, 1]) 



n2(n - 4)2 1 



16 



16 



(3.6) 



Let ^ — {\' -\- e)KQ, where e is arbitrary sufficient small, we need finally here 

16 ^^^^ 16 -p(A^ + g)^o>0- 
For that, it is sufficient to have for p — > +00 

- 8(e2 + e){n- 2)(n - 4) + o(-) > 0. 
16 p 

So (3.6) > holds only for n > 32 when p — )■ +00. Moreover, for p large enough 

8e2(n-2)(n-4) / + u^f^' < hJ < / |A<^r 

Thus it follows from Lemma 3.4 that u* is singular with A' = e^K^, /3 — {e^Ko + e{n,p)) 
and A* < e^Ko. 

(2) Assume 13 < n < 31. We shall show that u = ^3 5 satisfies the assumptions of 
Lemma 5.4 for each dimension 13 < n < 31. Using Maple, for each dimension 13 < n < 31 
one can verify that inequality (3.4) > holds for the A' given by Table 1. Then, by using 
Maple again, we show that there exists P > X' such that 

(n-2)2(n-4)2 1 



16 {\x\^-0.9\x\^+'){\x\^-\x\f 
^(n-l)(n-4)2 1 



> p^(l + W3.5)^+^ 



4 |xp(|xp — |x| 2 ) 

The above inequality and and improved Hardy-Rellich inequality (5.0) guarantee that the 
stability condition (5.2) holds for /3 > A'. Hence by Lemma 3.4 the extremal solution is 
singular for 13 < n < 31 the value of A' and /3 are shown in Table 1. 



Remark 1 The values of X' and (3 in Table 1 are not optimal. 

Remark 2 The improved Hardy-Rellich inequality (3.0) is crucial to prove that u* is 
singular in dimensions n > 13. Indeed by the classical Hardy-Rellich inequality and 
u := W2, Lemma 5.4 only implies that u* is singular n dimensions n > 32. 
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Table 1 



n 


A' 




31 


3MKo 


4.05Ko 


30-19 


4.6Ko 


lOKo 


18 


3.5Ko 




17 


3.26Xo 


3.60^0 


16 


3. 13^0 


3.78Xo 


15 


2.7QKo 


3.127^0 


14 


2MKo 


2.96/^0 


13 


2.03iro 


2.lhKQ 



